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We consider the gravitational model with additional spatial dimensions and anisotropic pressure 
which is nonzero only in these dimensions. Cosmological solutions in this model include accelerated 
expansion of the Universe at late age of its evolution and dynamical compactification of extra 
dimensions. This model describes observational data for Type Ia supernovae on the level or better 
than the ACDM model. We analyze two equations of state resulting in different predictions for 
further evolution, but in both variants the acceleration epoch is finite. 
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I. INTRODUCTION 

The most important event of last 15 years in astro- 
physics is conclusion about accelerated expansion of our 
Universe at late stage of its evolution. This conclusion 
was based on observations of luminosity distances and 
redshifts for the Typ_e Ia supernovae [J 0] , cosmic mi- 



large-scale galaxy clustering [1] , 



crowave background_ [3 1 
and other evidence 

To explain accelerated evolution of the Universe vari- 
ous mechanisms have been suggested, including the most 
popular cosmological model ACDM with a A term (dark 
energy) and cold dark matter (see reviews [HI])- The 
ACDM model with 4% fraction of visible (baryonic) mat- 
ter nowadays, 23% fraction of dark matter and 73% frac- 
tion of dark energy j.'jj describes Type Ia supernovae, 
data rather well and satisfies observational evidence, con- 
nected with rotational curves of galaxies, galaxy clusters 
and anisotropics of cosmic microwave background. How- 
ever, the ACDM model (along with vague nature of dark 
matter and energy) has some problems with fine tuning of 
the observed value of A, which is many orders of magni- 
tude smaller than expected vacuum energy density, and 
with different time dependence of dark energy SIa and 
material Q m fractions (we have Q\ ~ VL m nowadays). 

Therefore a large number of alternative cosmological 
models have been proposed. They include matter with 
nontrivial equations of state, for example, Chaplygin gas 
and modified equations of state fij-flil] ; scalar fields with a 
potential |12Hl4l |; modified gravity with f(R) Lagrangian 
[l5l [III; theories with extra dimensions, in particular, 
brane models [TtJ and many others d-Q. 

In this paper we explore the cosmological model with 
anisotropic pressure and nontrivial equation of state in 
1+3+d dimensions, suggested by Pahwa, Choudhury and 
Seshadri in Ref. jsj. The authors omitted the important 
case d = 1, we include it into consideration. We also 
analyse how to modify the equation of state and to avoid 
"the end of the world" (the finite-time future singularity) 



which is inevitable in the model [9|. 

In this model the 1 + 3 + d dimensional spacetime is 
symmetric and isotropic in two subspaces: in 3 usual 
spatial dimensions and in d extra dimensions. It has the 
following metric with two Robertson-Walker terms: 



ds 2 = -dt 2 



\t) 



dR 2 



dr 2 



1 — k\r 2 



r 2 dfl 



k 2 R 2 



R 2 dfl 



d-l 



(1.1) 



Here the signature is (—,+,...,+), the speed of light 
c = 1, a(t) and k\ is are the scale factor and curvature 
sign in usual dimensions, b(t) and k 2 are corresponding 
values for extra dimensions. It is supposed that the scale 
factor a(t) grows while b(t) diminishes, in other words, 
some form of dynamical compactification takes place. 

Authors of Ref. [§] develop the approach of Ref. [l8[ 
and suppose that the spacetime (|1.1|) is filled with a uni- 
form density matter with anisotropic pressure and the 
following energy-momentum tensor: 



T» = diag (-p, P a ,P a ,P a ,P b ,..., P b ). 



(1.2) 



Here p is the energy density and P a (Pb) is the pressure 
in normal (extra) dimensions. So in normal dimensions 
pressure is different from that in additional dimensions, 
while being isotropic within each subspace. 

In Ref. @ matter is supposed to behave like pressure- 
less dust (P a — 0) in usual dimensions, while in extra 
dimensions it has appreciable pressure Pb depending on 
density p by a power law 



(1.3) 



P b = Wp 1 - 



with a negative constant W. In this model matter (jTT 
with anisotropic pressure plays a role of dark energy and 
source of accelerated expansion. So the following Ein- 
stein equation without usual A term is considered: 



Gt = 8nGT,t. 



(1.4) 
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To describe accelerated expansion of the Universe au- 
thors of Ref. Q used cosmological solutions of their 
model with d > 2 additional dimensions. In Sect. ITTl of 



2 



this paper we show, that the case d — 1 is also applicable 
to describe acceleration. In Sect. lIIII we apply this model 
for all d > 1 to describing observational data for Type 
la supernovae and determine optimal model parameters. 
In Sect. |lV]we modify the model [t| to solve the above 
mentioned problem of "the end of the world" . 



II. COSMOLOGICAL SOLUTIONS 



For the considered metric (jl.ljl in the case = the 



Einstein tensor components G£ (p, v = 0,1, 
1 < i < 3 < 7) are [9): 



3, 



a b 



G o = _ 3d ^_ 3 !i._ 
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b 2 a 2 ' 
d(d - 1) b 2 



ki 



G\ = (l-d) 
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If we substitute these expressions into Eq. (|1.4|) and 
add the continuity condition T£„ = we obtain the sys- 
tem of cosmological equations. This system has the form 



a 

a b 
a b 
da 



a b 
a b 

■2 
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8ttG 

o, 

8ttG 



-Ph 



dt ipab) + Pba dt b 



= 0. 



(2.1) 
(2.2) 
(2.3) 
(2.4) 



in the case with d = 1 extra spatial dimension, that did 
not considered in Ref. @. Here pressure P a in "usual" 
dimension equals zero, as mentioned above. Eq. (|2.4I) is 
the continuity condition for d = 1 and P a — 0. 

Using the Hubble constant H ~ 2.28 • 10" 18 c _1 0] 
and the critical density 

Pc " 8^G 

at the present time, we make the following substitutions 



(2.5) 



r = Hot, 



P Pb A 1 a 

P=—, Pb = —, A = \og— , 



Pc 



Pc 



a 



B = \og- 
bo 
(2.6) 

and introduce dimensionless time r, density p, pressure 
pb and logarithms A, B of the scale factors (here an, bo 
are present time values of a and b) . 

We denote derivative with respect to r as primes and 
rewrite the system (|2.ip - (|2.4p as follows: 



2.-1 



A' + A'-B' - n fe e 
2A" + 3A' 2 + B" + B' 2 

A" + 2A' 2 - n k e 
p' + ZpA' + (p + p b )B' 



2A'B' = Q k e 

-2A 

" ~Pb, 

0. 



-2A 



(2.7) 
, (2.8) 

(2.9) 
(2.10) 



Here 



If we express 



n, 



-fci(ao-ffo) 



B' = {p + n k e- 2A )/A' - A' 



(2.11) 



(2.12) 



from Eq. (|2.7p and substitute it into three equations 
(HU) - (|2~T0)) . one should note that Eq. (gl§ may be 
reduced to Eq. (|2.9|) . So in the planar case 



ki=k 2 = 0, 



O fe = 



we have the system of two independent equations 

A" 
p' 



2A'-p b , (2.13) 
3pA' + (p+p b )(A' -p/A'). (2.14) 



If we fix an equation of state for pressure p b , for exam- 
ple, the above mentioned power law (|1.3[) 



Pb = wp 



1-7 



(2.15) 



we may consider the equations (|2.13[) . (|2.14[) as a closed 
system of first order differential equations with respect 
to 2 unknown functions A'(r) and p(r). The dependence 
(|2.15p is used in Ref. @, where parameters w and 7 are 
chosen in accordance with observations. 

The Cauchy problem for the system (|2TT31 . (I2TT4]) 
requires two initial conditions. We refer them to the 
present epoch (here and below it corresponds to the value 
t = 1) in the following form: 



A'l ,=1, 

I T— 1 5 



P\r=l = 



(2.16) 



The first condition results from definition of the Hubble 
constant 



a 



HaA'l 



In the second condition (|2.16[) we suppose that the energy 
density p = p ■ p c at the present time has the fraction fio 
in the critical density (|2.5[) . In Ref. [9] this fraction equals 
matter density fraction in the ACDM model 0: 



n = n m = 0.27. 



(2.17) 



Note that in Ref. [9| the second condition (|2.16p was used 
in the form p\ — 1, but the value ilo (|2.17p was taken 
as the factor in the r.h.s. of Eq. (jl.4l) . From our point 
of view, that approach introduces useless vagueness in 
physical sense of the value p. In our approach p in con- 
ditions (|2.16p is density of all gravitating matter (visible 
and dark) with described above anisotropic pressure. 

Remind that we have no dark energy or A term in 
Eq. (|1.4p in the model [9]. Anisotropic pressure in addi- 
tional dimensions plays here the role of dark energy as a 
source of acceleration. The contribution of this source is 
the term Qb = ~B'\ _. in the equality 



(2.18) 



3 



that results from equation (12.7[) , if we fix it at the present 
time r = 1. 

To obtain cosmological solutions for d = 1, k% = in 
this model we are to solve numerically the Cauchy prob- 
lem for the system (|2.13l) . (I2.14|) with initial conditions 
(|2.16|) moving into the past for t < 1 and into the fu- 
ture for r > 1. Then we integrate functions A'(r) and 
B'(t) (12.12)) keeping in mind Eqs. (|2.6|) and calculate de- 
pendence of the scale factors a = aoe A , b — b§e B and 
density p on dimensionless time r. 
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FIG. 1: Scale factors a,b, density p and acceleration parameter 
— q depending on dimensionless time r for Qo = 0.27, 7 = 0.9 
and specified values of w 

The results of calculation for scale factors a(r), 6(r), 
density p(r) and the acceleration parameter 



A" + A 



aa 

1 = = 7 



/2 



(2.19) 



(q is the deceleration parameter) are presented in Fig. [T] 
Here k\ = 0, fio = 0.27, 7 = 0.9 and 3 scenarios for 
w = —1.6 (dash-dotted line), w = —1.8 (solid lines) and 
w = — 2 (dashed lines) are shown. 

This evolution begins from infinite value of density p 
at some initial moment tq. We can see here two different 
variants for this beginning. For solutions with w = —1.6 
and w = —1.8 (we denote them as "regular" solutions) 
the scale factor a expands from a — like a ~ \J~t — tq 
at the initial stage whereas the scale factor b diminishes 



from initial infinite value up to values b ~ bo during some 
percent of total lifetime of this universe. This behavior 
of 6(r) looks like some variant of dynamical compactifi- 
cation, because the parameter bo is arbitrary one in this 
model, we may put bo to be sufficiently small. 

Another type of evolution ( "singular" solutions) is rep- 
resented with dashed lines in Fig. [T] for w — —2. For 
singular solutions infinite value of density p at t = To 
corresponds to nonzero value of the scale factor a and 
b = 0. Obviously, these solutions are nonphysical and 
should be excluded. 

All regular and singular solutions in Fig. [1] describe 
accelerated expansion (for the factor a) at late stage of 
evolution. Beginning of this stage may be seen in the 
graph of the acceleration parameter —q(r). Acceleration 
rate depends on the parameters w, 7, flo and the cur- 
vature fraction (12.111) Ofe = —ki(a Ho)~ 2 depending on 
the sign k\. If fife ^ (fci = ±1), one should use the 
system |J23j) - (|2T2l) instead of Eqs. (|2T3]l , ([2~T4|) . In 
this case we integrate numerically the function A'(t) si- 
multaneously with solving the Cauchy problem for the 
system (|2.9p - (|2.12p . We add here the natural initial 
condition -A| 1 = to conditions (|2.16p . 

For all reasonable values of four free parameters w, 
7, f2o, fife the stage of accelerated expansion appears to 
be finite, because density p inevitably vanishes in this 
model. In Fig. [1] this effect may be seen in the graphs 
p(r) with logarithmic scale in Y-direction. We denote 
the moment of zero density by t*: p(r„) = 0. For r > T. t 
density p becomes negative and nonphysical, all energy 
conditions (in particular, the weak energy condition) are 
violated. 

This finite-time future singularity may be classified as 
the Type IV singularity in accordance with the scheme 
from Refs. 0, Eaf. For this singularity a(r«) is nonzero, 
p(t*) equals zero, the effective density and pressure 



Peff 



3H 2 
8^G 



PcA'\ 



Peff 



2H + 3H 2 _2q-l 
8^G ~ 



■ - 



remain nonzero, but higher derivatives of H diverge at 

Note that the main features of the considered cosmo- 
logical solutions, in particular, the future singularity, fi- 
nite lifetime tq < t < t* and negative density for r > T. t 
take place not only for d = 1, but also for higher dimen- 
sions d > 2. In the case of d > 2 additional dimensions 
after substituting the components into Einstein equa- 
tion p.4[) and substitutions (|2.6I) in these equations and 
Eq. (|2.4[) we have in the flat case k\ = = the follow- 
ing system @, generalizing Eqs. (ESJj) - (12TT41 : 



A" = 
P 1 = 

B' = 



d(d - 1) B\\B' - A') - 3(d + 1) A' 2 - 3dp b 



d + 2 

-3pA' — d(p + pb) B', 



(2.20) 



y/3[(d + 2) A' 2 + 2{d-l)p]/d- 3 A' 



d-1 
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Solutions of the system (12.201) for d > 2 were obtained 
m Ref. @, but some features of them were not consid- 
ered in that paper. For example, singular solutions with 
nonzero value a(ro) (where p is infinite at the initial mo- 
ment To) also take place for d > 2, if the value w is less 
than the critical value w cr (l, ^o)- In Fig. [2] boundaries 
w = w cr separating domains of regular and singular so- 
lutions on the 7, w plane are presented for different d 
and £1 - Singular solutions are described by the inequal- 
ity w < ui cr (7, r2 ) and lie below corresponding lines in 

Fig.m 
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FIG. 2: Boundaries w — w cr between domains of regular 
(above) and singular solutions (below a curve) for indicated 
values Qo and d 

Another important property of these cosmological so- 
lutions is their finite-time future singularity in other 
words, inevitability of "the end of the world" because 
of vanishing density at r = t* for all d (see Fig. 0] below). 
The authors of Ref. @ did not pay attention to this phe- 
nomenon, essential for their model. It is connected with 
the chosen equation of state (|2.15[) for pressure pb in ex- 
tra dimensions. This drawback will be eliminated with 
modifying the model @ in Sect. IIVI after application this 
model to describing observational data for Type la su- 
pernovae in the next section. 

III. APPLICATION TO SUPERNOVAE 
OBSERVATIONS 

To apply the model to describing the observational 
data it is convenient, following the authors of [9(, to 
use Internet table [19j for Type la supernovae in distant 
galaxies. At the present moment this updated table con- 
tains redshifts z = Zi, distance moduli fii and errors Oi 
of ^ for N = 580 supernovae. 

Redshift 

2 = 4°, - 1=e -Mr)_ 1 (3.!) 

a[t) 

is associated with the value of a at the time t of a su- 
pernova light emission. The distance modulus fi is the 
logarithmic function 

10 pc 



of the luminosity distance @, 0] : 

o 1 
„ , , f dz On f df , 

Dl= 1 + 2 ) -57- = tt-TT /-FV 3 - 2 
J H(z) H a(T) J a[T) 

z r 

To describe the data [19| of Type la supernovae, for 
given values d, w, 7, Q,q of this model we consider evo- 
lution of the scale factor a(r) and dependence of the nu- 
merical integral (|3.2[) Dl and /1 on r. For each value of 
redshift Zi in the table [l9| we calculate the correspond- 
ing r = Tj with using Eq. (13.11) and linear approximation 
and the theoretical value fith = M 7 *) f° r r « fr° m Eq- (13.21) . 
The measure of differences between these theoretical val- 
ues [ith = Hth{d, w, 7, fioj ^fe, Zi) and the measured values 
Hi is Q: 

2 / j o o \ ~ Vthid, ■ ■ . ,Zi)] 
X [d, w, 7, S2 , S2fc) = 2^ 2 ■ ( 3 - 3 ) 

i=i 0i 

The authors of Ref. Q calculated optimal parameters 
w and 7, minimizing the function (|3.3I) for the flat model 
(fci = 0) with fixed ft = 0.27 (POT)) and d > 2. In this 
approach for each d > 2 they minimized the function 
X 2 (i«, 7) of two variables. 

We generalize their approach to the case d = 1 ad- 
ditional dimension. At the first step we fix k\ = 0, 
f^o = 0.27 in according with Ref. Q and obtain the pic- 
ture of level lines for the function \ 2 (, w il)i presented in 
Fig. [3] for d = 1 and d = 2. 



d=1 rf=2 




FIG. 3: Level lines of x 2 (™,7) for k% = 0, fi = 0.27. The 
dashed line is the boundary of singular solutions 

Here the dashed line is taken from Fig. [2] and separates 
regular and singular solutions. We see that for d — 1 and 
d = 2 the minimum of % 2 lies above this line, that is in 
the domain of regular solutions. The same picture also 
takes place for d > 3. 

For each d > 1 we calculated minimums for the func- 
tion of two variables x 2 { w il) an d coordinates w, 7 of 
this minimum. They are represented in Tabic |U 

We compare these minimal values with the value of 
the function (HP)) for the flat ACDM model with the 
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TABLE I: Minimum of \ 2 an( i optimal values w, 7 for fixed 
fci — 0, flo — 0.27 and various d 



d 


1 


2 


3 


6 


10 


ACDM 


2 

mm x 


563.15 


563.41 


563.52 


563.65 


563.7 


563.29 


W 


-1.77 


-1.355 


-1.212 


-1.067 


-1.01 




7 


0.925 


0.819 


0.772 


0.715 


0.686 





same parameters k\ — 0, TO = f2o = 0.27 (therefore, 
f^A = 0.73) and the same supernova data [l9|. We see 
that the predictions are rather close, and for d = 1 the 
model @ fits the data better than the flat ACDM model. 

At the next step for more precise estimation of opti- 
mal model parameters we consider variations of fractions 
fio = f^m and Jlfc for matter density and curvature re- 
spectively. One should take into account these degrees 
of freedom in both models: the model [§] and ACDM. 
In the model [9( for each d > 1 we minimize the func- 
tion (|3.3|) of four variables: x 2 (u>, 7, ^o> ^fe)- We also 
compare this results with the same value of the ACDM 
model (where \ 2 depends on Q m and Qk) and keep in 
mind the constraints on these parameters due to cosmic 
microwave background anisotropy, galaxy clustering and 
other factors Q: 



0.2743 ± 0.0072, -0.0133 < n k < 0.0084. (3.4) 



Numerical search of this minimum includes a start- 
ing point (for example, the values from Table HJ, analy- 
sis of gradients or increments for x 2 and the constraints 
(I3.4[) . The results of calculation with optimal values of 
the model parameters are presented in Table HU 



TABLE II: Minimum of x an d optimal values w, 7, Qo, &k 



d 


1 


2 


3 


6 


10 


ACDM 


min x 


563.136 


563.39 


563.506 


563.634 


563.69 


563.058 


w 


-1.740 


-1.323 


-1.2032 


-1.061 


-1.003 




7 


0.926 


0.821 


0.7739 


0.7174 


0.696 




n 


0.2815 


0.279 


0.274 


0.2673 


0.267 


0.2716 


fife 


0.0084 


0.0084 


0.0084 


0.0084 


0.0084 


-0.0133 



We see that the ACDM model is more sensitive to vari- 
ations of Q m and flk and the better result for this model 
is achieved. Here optimal values of the model parameters 
are determined by the constraints (|3.4[) . We impose these 
constraints on the model [|| though they are not strictly 
applicable to it. In this model minx 2 weakly depends on 
Qo and flk, so we can not diminish x 2 appreciably if we 



slightly broaden the limitations (|3.4j) . 

In Fig. |4] one can see evolution of the scale factor a(r), 
(and b for the model [9j), the acceleration parameter 
— q(r) and density p(r) for the ACDM model and the 
model [9( with d = 1 (solid lines), d = 2 (dots) and d = 6 
(dash-dotted lines). For all these models we use the op- 
timal parameters from Table [II] 
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FIG. 4: Scale factors ce(t), 6(t), acceleration parameter —q(j) 
and density p(r) for the optimal solutions from Table ITT1 



Evolution of the scale factor a(r) for the model Q 
with different d and for the ACDM model is very close 
up to z ~ 1.5 (a > 0.4ao), before this epoch the ACDM 
model demonstrates slower expansion. This difference is 
more visible for the acceleration graphs —q(r). The scale 
factor b for the case 0] diminishes to 6 ~ 6 according to 
the mentioned above compactification scheme (compare 
with the regular solutions in Fig. [l]). 

Behavior of cosmological solutions in the future for 
both models is also different. The ACDM model demon- 
strates unlimited accelerated expansion whereas for the 
model 0] the acceleration turns into deceleration and in- 
evitability results in the above mentioned zero density 
p at t = T,f with nonphysical values p < for t > t„. 
The finite lifetime of this universe depends on d, it is the 
smallest for d = 1. In the next section we discuss how to 
eliminate this essential drawback of the model. 



IV. MODIFICATION OF THE MODEL 

We have noted that all cosmological solutions in the 
model ^9f have the finite-time future singularity. This 
inevitable "end of the world" is connected with the cho- 
sen power law dependence (|2.15[) of pressure pb in extra 
dimensions on density p. The terms with the factor pb 
in equations (|2.14j) for d = 1 or (|2.20p for d > 1 deter- 
mine rate of density decreasing when p is small at the 
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end of its evolution. In this case the leading terms in the 
mentioned equations are 




p-> 0. 



(4.1) 



For p — > we have nonzero values A' and B', so for 
the weak power law dependence (|2 . 15[) the approximate 



equation (|4.1|) p 



-Cp 1 7 has the finite solution 

- [tC(-< 



I 1/7 



(4.2) 



To avoid this finiteness we are to modify the equation 
of state (power law dependence) (|2.15|) of the model @ 
for small p. In particular, a linear dependence for p close 
to zero 



Pb = w p, 



p^O 



(4.3) 



ensures infinite evolution with positive density. 

The linear law (|4.3|) for all p does not describe the ob- 
served accelerated expansion. For good agreement with 
observations we are to search an equation of state pb(p) 
with slower growth of \pb | at high p similar to Eq. (|2.15|) . 
We suggest the appropriate variant of this dependence 



Pb 



w 1 



Po 



(4.4) 



with the linear law (|4.3p for p <C po (here u>o = wi+w/po) 
and another linear law pb — wip for p ^> po- 

The model flU) - (|?T2"|) or f2~2"U| for d > 1 with the 
linear-fractional equation of state (I4.4[) makes it possible 
to avoid finite lifetime of the type (|4.2p and to transform 
it into the exponential asymptotic behavior 



p ~ exp(— Ct), C = const • ^u>i 



W 
Po 



(4.5) 



This behavior results from the equation p' ~ — Cp and 
may be observed in graphs p(r) in Fig. [5] 

For the model with Eq. (|4.4I) we can find optimal values 
of parameters w, wi, po, f^o, ^fc presented in Table Hill 
and achieve better agreement with the supernovae data 
19] than for the models ACDM and @ with Eq. (|2"T5l 
Cosmological solutions for the model with Eq. (|4.4[) and 
parameters from Table HTT1 are shown in Fig. [5] 



TABLE III: Optimal parameters for the model with Eq. (|4.4[1 
Po = 0.005 is fixed 



d 


i 


2 


3 


6 


10 


ACDM 


2 

mm x 


562.898 


562.814 


562.79 


562.766 


562.757 


563.058 


w 


-1.603 


-1.026 


-0.84 


-0.658 


-0.587 




Wl 


-0.195 


-0.343 


-0.387 


-0.426 


-0.439 






0.2815 


0.2716 


n k 


-0.0133 


-0.0133 



We see in Table IIIII that the accuracy of the model 
with Eq. (|4.4j) increases (x 2 diminishes) for large d, un- 
like in the case with Eq. (pT5l) in Table HI We should 
note that the values % 2 m Table IIIII are not absolutely 
minimal, because we fixed the parameter po = 0.005. It 
is interesting, that for all d we can achieve smaller values 
minx 2 , if we take smaller values of po. But if po — > 0, 



the factor C in the exponent (|4.5|) tends to infinity, the 
density p decreases too rapidly and the picture of van- 
ishing p looks like in the finite case in Fig. |4] So we put 
the restriction p > 0.005 to exclude this almost instan- 
taneous transition to the state with p ~ 0. Under this 
constraint we have the optimal value po = 0.005 and also 
flo = 0.2815, Q k = -0.0133 for all d. 

Fig. [5] demonstrates cosmological solutions for the 
model with Eq. (|4.4j) with the optimal values of parame- 
ters from Table UTTl For both models Eqs. and (|2~T5)) 
in Figs.[3]and[n]the acceleration epoch is finite and its du- 
ration depends on d in the same manner. But after this 
epoch for the model with Eq. (|4.4[) we see here infinite 
decelerated expansion. 
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FIG. 5: Scale factors o(t), b(r), functions — g(r) and p(r) for 
the model with Eq. (|4.4p with the parameters from Table Hill 



Graphs of a(r), — q(r) and p(r) here are more close to 
the dashed lines for the ACDM model during the accel- 
eration epoch than the similar curves in Fig. 2) But after 
the mentioned epoch predictions of the ACDM model and 
the model f)4.4[) sharply diverge. 

In Fig. [S] we present how the model with Eq. (I4.4[) for 
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d = 1,2,6 and the ACDM model describe the supernovae 
data from the site |19] (dots) in the z-Dl plane. All 
these models were considered below in Fig. [5j we use the 
same optimal parameters from Table IIIII and the same 
notations for the curves, in particular, the dashed line 
for the ACDM model. 




z 



FIG. 6: Luminosity distance Dl in Gpc depending on red- 
shift 2 for the models ACDM and (|4.4|l with parameters from 
Table HII] Dots are the data [jj] 

The presented curves are very close in the region z < 1 , 
for larger z the ACDM line slightly diverges from others. 
These lines are result of optimal fitting to the observa- 
tional data [H (580 dots in Fig. |BJ). The values x 2 in 
Table IIIII show rather good results for the model @ with 
Eq. (|4.4j) for pressure, but these values are not the best 
fit, because we fixed po to avoid the mentioned above 
sharp transition to p ~ 0. 

V. CONCLUSION 

The gravitational model @ with additional spatial di- 
mensions and anisotropic pressure provides accelerated 
expansion of the Universe corresponding to observational 
data for Type la supernovae . The authors of Ref. [1] 
did not consider the case with d = 1 additional dimen- 



sion, but we found that for the chosen in [9| power law 
equation of state (|2.15p it is the case d = 1 that yields 
the best fit (see Tables Hand HI. 

Unfortunately, the model |9[ with the power law de- 
pendence (|2.15|) inevitably predicts predicts the finite- 
time future singularity of Type IV in classification from 
Refs. @, EH . It is connected with vanishing density p at 
finite time r = (and negative density for r > r*). Evo- 
lution of this universe is broken at r* , the finite lifetime 
is shorter for small d (see Fig. |4j. 

We demonstrated in Sect. IIVI that this drawback has 
technical nature. It is connected with too weak depen- 
dence pb(p) for small p in the power law equation of state 
(|2.15p . If we modify this law and choose a linear de- 
pendence (I4.3P for small p, we obtain an infinite cosmo- 
logical evolution with positive density (but vanishing at 
t — y oo) . We suggest the linear- fractional variant (|4.4[) 
of dependence pb{p) to solve the following two problems: 

(a) to avoid "the end of the world" of the type (|4.2|) and 

(b) to describe 580 Type la supernovae data points from 
the site [l9| . The dependence (|4.4[) is a bit more compli- 
cated than Eq. (|2.15[) . but it successfully conserves pos- 
itive density p during infinite lifetime and fits the data 

19] better than the ACDM model and the model [|| with 
Eq. (J2T5]). 

Note that for the model with Eq. (14.41) we practically 
used only two parameters w and w± for fitting. The value 
Po was fixed because for very small po we have better fit, 
but the sharp downfall of p(r) to p ~ looks similar to 
the mentioned "end of the world" . The parameters £7o 
and Ofc influence on minimum of x 2 rather weakly for the 
model with Eq. (14.41) . If we fix, for example, Qq = 0.27 
and flk = 0, we obtain minimums for x 2 differing from 
results in Table iLUl less then 0.01 for all d. 

Cosmological solutions in the model with Eq. (14.41) are 
divided into regular and singular ones similarly to solu- 
tions with Eq. (|2.15l) shown in Fig. [T] However, Fig. [5] 
demonstrates that for the optimal values of parameters 
from Table 1X111 solutions with Eq. (|4.4[) are regular. 

It is interesting that the model [9j with both consid- 
ered variants of dependence pb on p (|2.15|) and (j4.4j) pre- 
dicts finiteness of the acceleration epoch. Its duration 
depends on d in the same manner (compare Figs. U]and 
[5]) and then acceleration sharply turns to deceleration. 
In the case (j2.15[) this evolution is broken at r = t* with 
/o(t«) = 0, but for the model with Eq. (j4.4l) the deceler- 
ated expansion is infinite and density p(r) tends to zero 
in the exponential form (|4.5[) . 
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